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Abstract. Directed paths have been used extensively in the scientific literature 
as a model of a linear polymer. Such paths models in particular the conformational 
entropy of a linear polymer and the effects it has on the free energy. These directed 
models are simplified versions of the self-avoiding walk, but they do nevertheless 
give insight into the phase behaviour of a polymer, and also serve as a tool to 
study the effects of conformational degrees of freedom in the behaviour of a linear 
polymer. In this paper we examine a directed path model of a linear polymer in a 
confining geometry (a wedge). The main focus of our attention is c n , the number 
of directed lattice paths of length n steps which takes steps in the North-East and 
South-East directions and which is confined to the wedge Y = ±X/p, where p is 
an integer. In this paper we examine the case p = 2 in detail, and we determine 
the generating function using the iterated kernel method. We also examine the 
asymtotics of c n . In particular, we show that 

c n = [0.67874...] x 2"- 1 (l + (-l)") + 0((4/3 3 / 4 ) n+ °( n )) + o ((4/3 3 / 4 ) n ) 

where we can determine the constant 0.67874 ... to arbitrary accuracy with little 
effort. 
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1. Introduction 

Lattice paths and lattice walks have long been used as models of conformational 
entropy in linear polymers O [10] . Perhaps the most famous of these models is the 
self-avoiding walk |221 116) . The most fundamental quantity in this model is w n , the 
number of self-avoiding walks of length n steps from the origin in the hypercubic 
lattice. Is it known that w n = pb n+ °^ n \ where /i is the growth constant of the self- 
avoiding walk, and that the limit linin^oo wl/ n — \i exists |12) . The self-avoiding walk 
is a non-Markovian model, and it is generally very difficult to extract its properties 
by rigorous or even by numerical means. 

Directed paths, in particular models in two dimensions, are generally recurrent 
models which can often be solved exactly by a renewal type argument. A particular 
example is the Dyck path which renews itself each time its visits the line Y = X in 
the square lattice. This property introduces a translational invariance in the model, 
which can be used directly to solve for the generating function as a root of a quadratic 
polynomial, see for example reference [8]. This general observation holds for other 
models, including a model of directed paths above the line Y = 2X [TT] and partially 
directed lattice paths including bargraph paths [21 [H [551 12H1 [2Z] and other related 
models such as Motzkin paths [5]. 

Confining a lattice path to a wedge in the square lattice introduces complexities 
which oftens makes the model harder to solve. These models are directed versions of 
the self-avoiding walk confined to a wedge [13] , which in turn is a model of a linear 
polymer in a confined geometry. A Dyck path is perhaps the simplest directed path 
model of a path confined in a wedge. In figure [lja) a more generic directed path in a 
wedge formed by the F-axis and the line Y — rX is illustrated (where r > 0). This 
model can be solved exactly for r £ {0, 1,2,3}, but no explicit solutions are known 
for other values of r although the radius of convergence of the generating function 
is known explicity for all r > 0, see for example references [71 [HI [HI [Ml [HI US]- 
Generally these models pose challenging combinatorial questions. If the line Y = rX 
has rational slope, then the generating function is a root of polynomial (it is algebraic), 
and a recurrance can be determined by a renewal type argument: The path renews 
itself each time its visits the line Y = rX. More generally r is irrational, and in these 
models there is no translational invariance of the model along the line Y — rX and a 
recurrence relation for the number of paths seems out of the question. 




The model in figure []Ja) can be generalised by considering a partially directed 
path above the line Y = rX as illustrated in figure [ljb). A recurrence relation for the 
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Figure 2. A fully directed path in the square lattice, starting at the origin, and 
giving steps in the North-East and South-East directions. 



generating function of these model has been written down for some values of r [18J, 
but generally these models pose a significant mathematical problem. A third variant 
of these models is illustrated in figure [2(c). In this case a partially directed path is 
included in the wedge formed by the X-axis and the line Y = rX. A model of this 
type was proposed and solved in the case that r = I in reference [21] . 

In this paper we generalise models of fully directed paths from the origin to a 
model of lattice paths in a wedge. Consider the directed path in figure [2] in the square 
lattice which takes steps only in the North-East and South-East directions. The most 
fundamental quantity in this model is c n , the number of paths from the origin of 
length n steps. Obviously, in this model c n — 2". This path can be put into a wedge 
generally as illustrated in figure [3] This is a symmetric wedge (about the X-axis), 
and is an alternative and more challenging model compared to the cases illustrated in 
figure QJ a). 

We are primarily interested in a model of directed paths starting from the vertex 
(2,0) and confined to the wedge formed by the lines Y — ±X/2. We determine 
recurrences for the generating function using the kernel method [24 , and in particular 
a variant of this method developed in reference [21] and which we will call the iterated 
kernel method. 



The most fundamental quantity in our model is c„, the number of paths of length 
n steps and confined to the wedge formed by Y — ±X/2 and starting from the vertex 
with coordinates (2, 0). We will solve for the generating function of c„, and our solution 
will be an alternating series of compositions of a root of a quartic polynomial. This will 



Y 




Figure 3. A directed path from the vertex (2, 0) in a symmetric wedge formed 
by the lines Y = ±X/2. 
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allow us to determine c„ to high accuracy: in particular, we show that asymptotically, 

c n = [0.67874 . . .] x 2"- 1 (l + (-1)") + O ((4/3 3/4 ) n+o(ll) ) + o ((4/3 3/4 ) n ) (1) 

where the constant 0.67874 can be determined to hundreds of significant digits with 
minimal computational effort. 

In section 2 we define our models and determine a recurrence for the generating 

1 /n 

function. In addition, we note that lim^-joo c n — 2, we solve for the generating 
function using the kernel method, and examine the properties of the roots of the 
kernel. These results allows us to compute the constant in equation |T]) in section 3 by 
examining the singularities in the generating function. We conclude the paper with 
some final comments in section 4. 

2. Directed Paths in a Wedge 

Let Ej be the square lattice of points with integer coordinates in the plane. A directed 
path in this lattice is a path that takes only North-East (NE) and South-East (SE) 
steps. If the path consists of n steps (or edges), then there are 2" such paths. One 
such path is illustrated in figure [U Let XY be the usual Cartesian coordinate system 
in figure [2J with origin at the first vertex of the path. Then the edges in the directed 
path have length y/2 each. 

The directed path in figure [2] is unconstrained by the boundaries of the wedge 

Y = ±X. This model becomes more interesting if the path is constrained by a 
narrower wedge Y — ±X/p, where p > 1 is an integer. 

Define the 1/p- wedge V p by 

V p = {(x, y) G Ej I where — x/p < y < x/p}. (2) 

Then V p is the subset of TL in the first and fourth quadrants bounded by the lines 

Y — zbX/p. In figure [3] a directed path confined to the wedge V p is denoted where 
p = 2. This path has its first vertex at the point with coordinates (2,0). Generally, 
directed paths in the wedge Y — -bX/p will have their first vertices at (p, 0). 

2.1. Directed paths in the Wedge formed by Y — ±X/2 

In this section we determine a recurrance for the generating function G2 of directed 
paths giving NE and SE steps from the vertex (2,0) in the wedge V2. Proceed by 
introducing the edge generating variable t and define 

G 2 (a,b)= c ( V(u,v)u a v b t n (3) 

n,iz,?;>0 

where c n (u, v) is the number of directed paths in the wedge V2 from the vertex (2, 0) 
of length n and with final vertex a vertical distance u from Y = \X/2~\ and a vertical 
distance v from Y = — \X/2] . 

Observe that the unit of vertical distance is determined by the fact that the 
squares in figure [3] have diagonal length 2. For example, a NE-step will increase 
the F-coordinate of the endpoint by 1, and the distance between the endpoint and 

Y = [X/2] may either decrease by 1 if the initial X-coordinate is even, or remain 
unchanged otherwise. 

Consider the examples in figure 2] The vertical distances are measured to the 
step- functions Y — ±\X/2~], and in these examples the vertical distances are u = 11 
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and v = 7 (figure HJa)) and u — 13 and v = 7 (figure HJb)). In addition, u + v is 
always even. 

(2) (2) 

We are fundamentally interested in c„ = ^2 u v>Q c n (u, v), the number of paths 
of length n in the wedge V2. The generating variables a and b are introduced to 
account for the vertical distance of the endpoint from the wedge boundaries and to 
enable us to determine G2(l, 1)- 

The examples of paths in figure 2] illustrates that the enumeration of these paths 
are fundamentally affected by a parity effect. Paths of even length (figure Ufa)) are 
either the single vertex at the point X = 2 generated by ab, or can be generated by 
adding a NE-step or SE-step to a path of odd length. Paths of odd length can only 
be generated by extending a path of even length by the addition of a NE-step or a 
SE-step. 

Let go(a,b) generate paths of even length, and suppose that gi(a,b) generates 
paths of odd length. Then 

G 2 (a,b) = g (a,b) + gi(a,b). (4) 

Thus, by determining go(a, b) and g\ (a, b), G2(a, b) can be determined. 

Recurrences for go(a,b) and gi(a,b) can be obtained by arguing as in figures [5] 
and[5J The basic idea is to create paths generated by <7i(a, b) by appending an edge to 
paths of even length generated by go (a, 6), and to create paths generated by go(a,b) 
by appending an edge to paths generated by g\ (a, b). The resulting recurrences are 

g {a,b) = ab + t{a/b + b/a)gi{a,b) - t(a/b)gi(a,0) - t(b/a)gi(0,b), 

9l (a, b) = t(a 2 + b 2 )g {a, b) - ta 2 g (a, 0) - tb 2 g Q (0, b). (5) 

These recurrences can be iterated to enumerate the directed paths. The numbers for 
even length paths are given in table 12.11 

The number of digits in the c„ in table 12.11 increases linearly with n, which 
implies exponential growth of c„. It is in fact possible to prove explicitly that Cn 
increases proportional to 2", using the techniques in reference 13J . In addition, since 
CnCm < Cn+m (concatenate a path of length m with a path of length n by translating 
the first until its first vertex coincide with the last vertex of the second to see this), 
c„ is a supermultiplicative function of n, and thus the limit 

lim c\l n = 2 (6) 

n — *oo 

exists 14 . Thus, there is a function of n, Co = e°^ n \ such that 

C = C 2" + o(2"). (7) 
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t(b/a) t(a/b) 



9o(a,b) = 




Figure 5. Determining go(a,b). Each path generated by go(a,b) is either the 
single vertex at (2, 0) generated by ab, or it is generated by appending a NE- 
edge onto a path of odd length (this gives the term (tab)(b/a)gi(a,b)), or it 
is generated by appending a S-E-edge onto a path of odd length (giving the 
term (tab)(a/b)gi(a,b)). Lastly, paths which step outside the wedge must be 
subtracted: (tab)(b/a)gi (0, b) if the path steps over the line Y = X/2 and 
(tab)(a/b)gi (a, 0) if the path steps over the line Y = 

Examination of the data in table 12.11 shows that c„/2™ approaches a constant. 
Assuming that Co is a constant, one may estimate it numerically. Dividing c n by 
2™ and increasing n up to n = 82 gives Co to five digits, namely 

C = 0.67874 ... (8) 

We proceed next by computing Co by solving the recurrances in equations ([5|) 
above, and by examining the singularities in gv{a,b). 

2.2. Solving equations {5p 

In this section we use the iterated kernel method (see reference [21]) to find an 
expression for the generating function go(a, b). Simplify the recurrences by introducing 



(tab)(b/a) (tab)(a/b) 



9o(a,b) = 




o) 



Figure 6. Determining gi(a,b). Each path generated by gi(a,b) is obtained 
by appending a NE-edge or SB-edge onto a path of even length. In this case, 
the vertical distance to the wedge does not increase by a full step, and thus no 
new factors of a or b are included, apart from accounting for stepping closer 
to the top boundary if a NE-edge is added, or stepping closer to the bottom 
boundary if a SB-edge is added. Thus it generates the term t(a/b+b/a)go(a, b) as 
indicated above. Lastly, paths which step outside the wedge must be subtracted: 
t(b/a)go(0,b) if the path steps over the line Y = X/2 and t(a/b)go(a,0) if the 
path steps over the line Y = —X/2. 
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n 




n 







l 


42 


2985401474160 


2 


4 


44 


11941093593120 


4 


12 


46 


47764374372480 


6 


48 


48 


191053247884320 


8 


180 


50 


764212991537280 


10 


720 


52 


3056816328436200 


12 


2820 


54 


12227265313744800 


14 


11280 


56 


48908759609676540 


16 


44760 


58 


195635038438706160 


18 


179040 


60 


782537580134560920 


20 


713760 


62 


3130150320538243680 


22 


2855040 


64 


12520579171583415840 


24 


11403060 


66 


50082316686333663360 


26 


45612240 


68 


200329075631136029040 


28 


182321460 


70 


801316302524544116160 


30 


729285840 


72 


3205263549296411867340 


32 


2916160800 


74 


12821054197185647469360 


34 


11664643200 


76 


51284202287042290859820 


36 


46650808680 


78 


205136809148169163439280 


38 


186603234720 


80 


820547109423871153955280 


40 


746350368540 


82 


3282188437695484615821120 



Table 1. The number of directed paths in V2. 



K = t(a/b + b/a) and L = t(a 2 + b 2 ). Observe also that go(a, 0) = go(0, a), and that 
9i (a, 0) = 51 (0, a). 

Since paths of odd length generated by gi(a,b) cannot intersect the boundaries 
of the wedge, they always will be weighted by factors of a and b. Hence, one expects 
gi(a, 0) = 5i(0, b) — identically. Ignoring this last observation for the moment gives 
the recurrances in slightly simplified form: 

go(a, b)=ab + K 9l (a, b) - t(a/b) gi (a, 0) - t(b/a) 9l (b, 0), (9) 

gi (a,b) = L g Q (a,b) - ta 2 go{ ai 0) ~ tb 2 g (b,0). (10) 

Substitute these equations into one another, and write them in kernel form. This gives 

(l-KL)g (a,b) 

= ab- ta 2 Kg {a, 0) - tb 2 K g (b, 0) - t(a/b) gi (a, 0) - t{b/a) 9l {b, 0), (11) 
(l-KL) 9l (a,b) 

= Lab - ta 2 g (a 1 0) - tb 2 g (b, 0) - t(a/b)L gi (a, 0) - t(b/a)L 31(6, 0). (12) 

We identify the kernel (1 — KL) in these equations. Generally, we say that recurrence 
equations are in kernelized form if the generating function and its coefficients have 
been collected on the left hand side, while all other terms and boundary terms are 
on the right hand side. The kernel (1 — KL) may be simplified, and then the quartic 
polynomial 

t 2 (a 2 +b 2 ) 2 -ab (13) 
in a and b appears as a factor. 
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To proceed, consider this to be a quartic in b with a and t two parameters. To 
solve the original recurrances, the roots of this quartic must be determined. Closer 
examination shows that the four roots have the following properties: The first is f3o( a ), 
which is a power-series in t: 

Pa(a) = a 3 t 2 + 2a 7 t 6 + 9a n i 10 + 52a 15 t u + ... (14) 

The second real root is /3i(a), which is singular at t = 0: 

„ , . a 1 / 3 2a b 'H 2 ' 3 28a 13 / 3 t 10 / 3 , s 

A(a) = ^ § si— + - (15) 

while the two remaining roots are a complex conjugate pair: 

R( . fl 1 / 3 (l T »V3) « 5 / 3 t 2 / 3 (l± ? V3) aH 2 

^( fl ) = ^3 + 3 3- + -'- (16) 

Substituting & = (3o(a) = (3q in the kernelized equations in ([TT|) and (fl2|) produces 
two equations for go(a, 0), <7o(&, 0), gi(a, 0) and gi(&, 0): 

K (ta 2 g (a,0) +t(3 2 g {(3 ,0)) + (t(a/(3 )gi{a,0) +t((3 /a) gi (l3 ,0)) = afo, (17) 

(ta 2 ffo (a,0)+^ 2 3 o(&,0)) +L(t(a/A))ffi(a,0)+i(A)/a)5i(A),0)) = Laf3 . (18) 
The solution of this linear system is 

ta 2 g {a, 0) + tf3 2 g (p , 0) = La/3 Q , (19) 

t(o/A))5i(a,0)+t(A)/a)ffi(A),0) = 0. (20) 
From the second of these it follows that 

gi (f3 (a),0) = -^ 9l (a,0) (21) 
Po 

and since this generating function cannot be negative, the conclusion is that 

3i(a,0)=0 (22) 
identically, as claimed above. The first solution above gives 

L( a,f3 )(3 _ ffi 
ta a 2 

Proceed by defining (a) = (/3o ° A) ° • • • A)) ( a ) to be the composition of (3o 
n-times with itself. Define (3^ (a) = a, then the last equation may be written as 

MP »«)- ^(n-i) I ^3iy I W ,0). (24) 

This may be iterated to obtain a solution for go(a,0) = go ((3^ , 0) : 

5oM) = i£(-iri(^ n) )^ n) /3^ +1) - (25) 

n=0 

Thus, one may solve for the generating function go(a, b) from the kernelized equations 
in IQTJ and J^J: 

o6-jg (a 2 go( a ,0) + 6 2 g (6,0)) 
0o(a,&) = — 1 _ RL (26) 

and the radius of convergence is given by the dominant root of the quartic in equation 
(|13p . provided that go(a,0) does not have compensating singularities at the same 
points. 



ffo(a,0 ) = ; j3o(/3o,0 . 23 
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2.3. More on the roots oft 2 (a 2 + b 2 ) 2 - ab. 
One may check that 

Ma) = [2ta]^\^ n ^ (27) 

n— v / 

as given in equation (fT4")l . The radius of convergence of this series is |at| 4 = 3 3 /4 4 . 

Closer examination also shows that fto(a) counts directed paths with first step 
in the SE-direction, above the line Y = —X/2 and with last vertex in the line 
Y = —X/2. The root (3i is the inverse function of [3 : Direct calculation shows 
that Pq o (3\{a) = j3\ o f3 (a) = a. 

The roots /?+ and /3_ are not independent, but one may check that 

c 3 \ „ ( —ix 3 \ ix it 2 / 3 x 5 it 2 x 9 7it w / 3 x 13 



^{-)=-P-{—)=-2l^ + — + T^ + ^6^ + --- (28) 

Solving for the roots of the quartic is equivalent to solving the non-linear system 

r 2 =(a 2 +b 2 ), (29) 

t 2 r 4 = ab (30) 

for a and b. There are four solutions, two given by the pairs (ai(r), b\(r)) and 
(— ai(r), — b\(r)), where 



, v V2r 3 t 2 ry/l -y/1- 4t 4 r 4 

aiW = , ; = = Ik 5 (31) 

VI + \/l - 4t 4 r 4 V2 

, , , r V / l + v / l-4i 4 r 4 v/2r 3 t 2 

V2 Vl - Vl - 4t 4 r 4 

and two more given by the pairs (&i(r), ffli(r)) and (— &i(r), — ffli(r)) (where we 
interchanged a\ and &i). An expression for r 2 is given below in equation (|35|) . One 
may check as well that 

oi (r) = tV f 1 + tV V f 4 " + 3 ) - J (33) 

u ; V froV 2n / (2n + i)4»+i/ 2 y ^ ; 

/ 4 4 /4n+ 1\ (tr) 4 " \ , . 

&i(r )=r l-t 4 r*> ( „ ,» ^, 34 



The roots of the quartic may be found by inverting ai to obtain r a (a) so that 
r a o aj = a\ o r a is the identity map. Then fto(a) = b\ o r a . Inverting b% to obtain 
rj,(6) gives a second root by the composition f3\ — a\ o r^. In particular, this means 
for example that (3q o fi\ = b\ o (r a oaijort = b\ o r& = identity since a^f 1 = r a and 
o^ 1 = rb- This proves the observation above that Po°Pi = Pi Po is the identity map. 
In other words, the composition of two roots of the quartic is the identity. 

The other two roots of the quartic is given by the compositions b\ o t\, and a\ o r a . 
Unfortunately, while explicit expressions for r a and rb can be obtained, they are rather 
lengthy; both r 2 and r\ are roots of the quartic t 4 x 4 — c 2 x + c 4 where c — a for r a 
and c = b for r&. This may be examined by iteration to determine the first few terms 
in r 2 . Comparison of the results to the online encyclopedia of integers [28] shows that 

2/ \ 2 f f4n\ (at) 4 " 

n— x 7 



Directed Paths in a Wedge 



10 



The second root of t A x A — c 2 x + c A proposes the " unphysical" series starting with 
(a/t 2 ) 2 / 3 + 0(a 2 ) for r 2 (a). The series expression for r 2 (a) may finally be substituted 
in &i(r) to obtain an expression for the root /3o( a ) of the quartic: 



A) (a) 



r a (a)y/l + y/T-4t*ri(a) 



V2 



(36) 



Remarkably, this evaluates to equation (|2T[) and compositions of this with itself will 
eventually lead to the expression for g Q (a,0) in equation (|2"5")) . 

In addition, having determined /3q (a) , one may consider the composition of a\(r) 
and r a (a), which is 



r a {a)^l-^T-At A r A a (a) 



V2 



(37) 



and which must be the identity map. In other words, by appealing to equation (|30p 
it follows that 



a ■ p (a) = t 2 r a (a) A , 
and from equation (|27p one concludes that the identity 

'4n\ (a<) 4 " 



E 

.71=0 



n / 3n + 1 



4n+l\ (at) 
n 



4)1 



3n + 2 



(38) 



(39) 



should be closely related to the combinatorial properties of directed paths in the wedge 
V 2 . 



3. Determining Co- 

In this section we examine the generating function go(a,b) in equation Ij26p with 
a = b — 1. Singularities in this generating function arise from several possible sources. 
In the first instance there are simple poles at the zeros of the kernel (1 — KL) in the 
denominator. These are located at t = ±1/2. In order to determine the constant Co, 
we examine the residue of go(l, l)/i n+1 at t = ±1/2. 

Put a = 1 and t = 1/2 in 0o(a). Compositions of 0o(a) with itself at this point 
gives the following values 





(1) 


= 1 








(1) 


= 2.9559774252208477098 . 


. x 10" 


-l 




(1) 


= 6.4633625443847777820 . 


. x 10" 


-3 




(1) 


= 6.7501832073150278963. 


. x 10" 


-8 




(1) 


= 7.6892979457392165146. 


. x 10" 


-23 




(1) 


= 1.1365801752937162161. 


. x 10" 


-67 




(1) 


= 3.6706268625677440729 . 


. x 10" 


-202 



where explicitly 

A)(l) = ^7^(19 + 3V33) 2 / 3 + ^^(19 + 3V33) 1 / 3 - \, 
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and the other terms are more complicated expressions involving nested radicals which 
do not simplify to manageable expressions. 

Generally, we observe that /3g (1) = a n at t = 1/2, and one may check that 
q„ + i w al/4. For example, a|/4 = 7.689297941 . . . x 10~ 23 w a 4 . This observation 
follows from the expansion /3o(a) = a 3 ^ 2 + 0(a 7 t 6 ) so that compositions of /3q at 
t = 1/2 and a = 1 quickly converges to zero. In other words, the recurrence 
x n+ i = (3o(x n )\t=i/2 is a fixed point iteration of order three. This fast convergence 
allows the accurate numerical estimation of /3q (1) at t = 1/2. 

At a = b = 1 the expression of go(l, 1) in terms of go(l,Q) and <?o(0, 1) is given 

by, 

^.^ ijMteMl , (40) 



where 



30(1,0) =30(0,1) = V(-l)"L(/3^)/3^/3(" +1 \ 



n=0 



and where L(a) = a 2 + /3 2 . Numerical evaluation of the residue at t = 1/2 using the 
calculated values of f3^ n '(l) above, gives the leading order behaviour of the number of 
paths of length n: 

c n = 2"- 1 x 0.6787405307981094574172327 . . . + parity term + . . . 

where the seven values of /3q (1) listed above produces C accurately up to at least 
O(10 -401 ) or 400 digits if each is calculated to at least this accuracy by using equations 
(|38p or (|27p or by explicitly using the closed form expression for /3q (a) and computing 
it to high accuracy using a symbolic computations package such as Maple. 

To determine the parity effects, we determine the residue at the pole located 
at t = —1/2 by repeating the analysis above. Put a = 1 and t — —1/2 in (3q(o). 
Compositions of fio(a) with itself at this point gives the identical values obtained 
above for t = 1/2. Thus, we conclude that 

c„ = 2 ,l ~ 1 (l + (-1)™) x (0.6787405307981094574172327. . .) + ... 

This is not an unexpected result since g$ (1, 1) should only enumerate even length 
paths; we note that C2 n +i = in the generating function go(l, 1). This result verifies 
the numerical estimate for Co from the data in table 12.11 in equation |(5J| . 

Contributions to c n also arise from singularities in the numerator in equation 
(f2Hj) . In particular, there are branch points and possibly other singularities in go(l, 0) 
and 3o(0, 1), and these are due to branch points in Pq. 

The radius of convergence of (3q can be determined when a = 1 by examining 
equations (|35| and (|38j) with a = 1. In particular, fi(l) in equation ([35| is convergent 
for all \t\ < 3 3 / 4 /4; in fact, evaluation shows directly that ri(l) is convergent for 
|t| = 3 3/4 /4. Since 3 3/4 /4 > 1/2, this proves that the simple poles at t = ±1/2 are 
within the radius of convergence of Po(l). 

Further examination of (3q (a) by Maple shows a complicated expression of nested 
radicals which expicitly contains factors of the form \/81 — 768a 4 £ 4 . This shows that 
there are branch points at at = 3 3 / 4 lj / 4 where u is a fourth root of unity. There may 
be more branchpoints on the circle at — 3 3 / 4 /4, but we did not verify this, and this 
will not play a crucial role in what follows. 
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Next, we consider the compositions of /3 in the alternating sum definition of 
<7o (1, 0). Since r\(f) is a positive term power series it follows by the triangle inequality 
in equation (|3"5)) that for (complex) t such that \t\ < 3 3 / 4 /4, 



Thus, by equation (j38|) 



(3 3 / 4 /4) 4 " _ 4 
3n + 1 ~~ 3' 



(41) 



\Po(l)\ < [3 3/4 /4] 2 (16/9) = if |i| < 3 3 / 4 /4. (42) 

Since (3o(a) is a power series with positive coefficients in both a and t, \Po(a)\ < |/3o(l)| 
for any \t\ < 3 3 / 4 /4 and \a\ < 1, and by the triangle inequality, /3o(a) is a maximum 
when t — 3 3 / 4 /4 for a fixed value of a. Thus, for fixed values of \a\ < 1, Po(a) is a 
maximum on the closed disk \t\ < 3 3 / 4 /4 when t = 3 3 / 4 /4. 

The above shows that for t = 3 3 / 4 /4, \0o(a)\ < l/y/3 for |o| < 1. Moreover, it 
follows from equation (gT)) that \0o(a)\ < M/a/3 for |o| < 1 when t = 3 3/4 /4. In other 
words, 

|A)(A)(1))I < |3 3 / 4 /4V3(V3 • 3 3 / 4 /4)| = -i. (43) 
It follows inductively that 

^ n) (l)|<^r. (44) 

Since the branch points in /3q"^(1) will occur at values of t that |/3q" (1)| = 3 3 / 4 /4, 
and < 1/V3" -1 < 3 3/4 /4 at a = 1 and for n > 2, this implies that the 

branch points in f3^ n \\) for n > 1 lies outside the circle with radius \t\ = 3 3 / 4 /4, and 
contributions of these branch points to c n are dominated by the contributions to c„ 
which is a result of the branch points in /3o(l) itself. 

In particular, since the radius of convergence of (3o(a) is determined by \at\ = 

3 3 / 4 /4 and |^ Tl_1) (l)| < l/VT" 1 , it follows that the radius of convergence of /^ n) (l) 

in the i-plane is on or outside the circle with radius \t\ — 3 3 / 4 v / 3"' 1 /4, for n > 2. 
Thus, the generating function go(a,0) has infinitely many singularities in the <-planc 
for a = 1, and by remark 1 following proposition 9 in reference [2j, go{l, 0) cannot be 
holonomic (or D-finite). Thus, go(l, 1), the generating function of even length directed 
paths in the wedge formed by Y = ±X/2, is not holonomic. 

The bound in equation ([44]) also proves that go(l, 0) is an absolutely convergent 
series on the open disk with radius \t\ — 3\/3 which includes the simple poles of go (1, 1) 
at t = ±1 and the branch points on the circle \t\ = 3 3 / 4 /4 in its interior. In other 
words, the asymptotic behaviour of c n is given by 

c n = 2"- 4 (l + (-1)") x (0.6787405307981094574172327 . . .) 

+ corrections. (45) 

The corrections are due to the branch points in /3q(1) and they grow at the 
exponential rate (4/33/4)n+o(n) _ Since 4 / 3 3/4 _ 1 ?5 < 2} the effects of the 
correction terms will dissappear fast with increasing n, and c n /2™ will approach 
0.6787405307981094574172327. . . at an exponential rate with increasing n. 
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4. Narrower Wedges 



It is possible to consider this problem in the narrower wedge V p where a and b measure 
vertical distance to the functions Y = ± \X/p] and where the path starts at the vertex 
with coordinates (p, 0). In this case the paths are counted by a generating function 
G p (a,b) = X)f=o 9i( a 'b) where gi(a,b) generates paths of length imodp and satisfies 
the set of coupled recurrance equations with a derivative term: 



,9o 

91 

92 
.93 



ab + t(a/b 
i(a 2 H 

t{a/b 



b/a) g p _ 



b 2 )g - ta 2 5o(a,0) 
- b/a) gi — ta 



tb 2 g (0,b), 



t(a/b + b/a) g 2 — ta 



~dgi 




-tb 


'dg{ 


db _ 


b=0 


da 


~dg 2 ~ 




- tb 


~dg 2 


db _ 


6=0 


da 



a=0 



~dg p -2~ 


-tb 


~dg p - 2 ~ 




db 


6=0 


da 


a=0 



g p -i = t(a/b + b/a) g p - 2 - ta 

Putting p = 2 in the above recovers the recurrences in equation ([5]) , if the fact that 
gi (a, 0) =51(0,6) = is used. 

Putting a = b = 1, writing G = J^f = g gi an d defining the boundary terms 



6=0 



for j = 2,3,. 



1 



Ag o (a,0) = 0, A 9l {a,0) = g x (a, 0) and A gj (a,0) 

and similarly for Agj(0, b) 7 show that by summing the above recurrances one obtains 



9gj 



p-1 



G=l + 2tG-J2 (&9i(h 0) + A 3 ,(0, 1)) . 



(46) 



The boundary terms J2i=o (A.^^ (1, 0) + A^.^O, 1)) represent the generating function 
of paths starting at the point (p, 0) and ending within one step from the lines 
Y = ±X/p. Such paths of length 2pn must give almost (within a constant) 
(p + l)n North East steps. The number of these paths grows at the exponential 
rate A = ((2p) 2p /(p— l) p ~ 1 (p+ l) p+1 )™. Hence the generating function of these paths 
is convergent inside the circle of radius l/A 1 / 2p . One can show that this is strictly 
greater than 1/2, hence all the boundary generating functions A<7i(l,0) and Ag^O, 1) 
of paths ending near the boundary of the wedge are convergent on the disk with radius 
strictly greater than 1/2. 

Thus, rearranging equation ([4^]) to get 

r l-E£o 1 (Ag,(l,0) + A. 9t (0,l)) 

G - T^2t (47) 

the dominant singularity in G is at t = 1/2 and is a simple pole. The subdominant 
singularities give terms growing to exponential order A" +0 *- n - ) where 



In other words. 



A = ((2p/(p - \) { 'P- l ^ 2 'P{p + i)(p+i)/2p)«) 
c ( n p) = A p 2 n + 0(X n ). 



(48) 



(49) 

Putting p = 2 produces the leading order term and the order of the next term obtained 
in equation ([45]). 
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5. Conclusions 

The main results in this paper are given by equations (|26|) and (|45j) . We have shown 
that the model of fully directed paths in the wedge V2 in figure [3] can be solved using 
the iterated kernel method. This model is the second in a sequence of models in the 
wedge V p , for p € {1, 2,3,.. .} and recurrences for p > 2 is given in section 2.4 (when 
p = 1 the model is trivial, and c„ = 2 n ). We proved that the number of fully directed 
path of length n in V 2 increases exponentially at the rate C$ 2™, plus subdominant 
terms which are also exponentials. In addition, we showed that Co (the coefficient of 
the leading order term) can be determined to high accuracy. 

In the more general wedge V p , it may be shown, using the techniques in reference 
[13], that the number of paths of length n is = Cq 2 n + lesser terms. In this event 
C^ p) is known for p = 1 (C { 1] = 1) and to high accuracy forp = 2 (C { 2) = 0.6 7 8 7 40 . . .). 
No estimates exist for other values of p. 

The "physical" root (3o( a ) of the quartic kernel (see equation (|27[) ) counts directed 
paths from (2, 0) (see figure [3]) above the line Y = —X/2 and with last vertex in 
the line Y = —X/2. By reversing the horizontal direction, this is also equal to the 
number of paths from the origin above the line Y — X/2 and with last vertex in the 
line Y = X/2. The generating function in equation (f25|) is an alternating series of 
products of compositions of 0o(a), suggesting that that go(a, 0) is "constructed" by an 
inclusion-exclusion process of directed paths above the lines Y = —X/2 and Y — X/2. 

These observations suggest that a combinatorial explanation for the generating 
function go(a,b) in equation (|26p might be possible in terms of paths "bouncing off" 
the boundaries of the wedge V2. We have not been able to find such an explanation. 
We leave this question open, but do observe that an explicit explanation in this case 
might give clues about the nature of the generating function of models in more general 
wedges, including the models of directed paths in the wedges V p for p > 2 formed by 
the lines Y — ±X/p, as well as in the half-wedges formed by the X-axis and the 
lines Y = X/p. Such models in half-wedges will be more difficult to solve, since 
the symmetric nature of V% plays a key role in writing down equation (|26|) and then 
determining Co by examining the singularities in equation (|40|) . 

Finally, our models in the wedge V p map to the well-known problem of random 
walks with given step-sets in the quarter plane (the first quadrant). These have been 
studied for example in references [TJ [2J [23] . The models in this paper correspond to 
models with starting vertex (p — l,p— 1) and step set {(1 — p,p + 1), (p + 1, 1 — p)}. 
For p = 2, this produces the generalised knight's walk model with step set (—1,3) 
and (3,-1) [21]. The knight's walk problem itself (with step set {(-1, 2), (2, -1}) 
and starting vertex (1,1) was studied in reference [2], where it was shown that its 
generating function is not holonomic (or D-finite). 

Acknowledgments: EJJvR and AR are supported by Discovery Grants from 
NSERC (Canada). 

Bibliography 

[1] Bousquet-Melou M 2001 Counting Walks in the Quarter Plane. In Mathematics and Computer 
Science: Algorithms, Trees, Combinatorics and probabilities, Trends in Mathematics 49-67 
(Birkhauser). 



Directed Paths in a Wedge 



15 



[2] Bousquct-Mclou M and Pctkovsck M 2003 Walks Confined in a Quadrant are not always D- 

finite. Theoretical Computer Science 307(2) 257-273. 
[3] Brak R, Guttmann A J and Whittington S G 1991 On the Behaviour of Collapsing Linear and 

Branch Polymers. J. Math. Chem. 8 255-267. 
[4] Brak R, Guttmann A J and Whittington S G 1992 A Collapse Transition in a Directed Walk 

Model. J. Phys. A: Math. Gen. 25 2437-2446. 
[5] Delest M-P and Viennot X G 1984 Algebraic Languages and Polyominoe Enumeration. Thcor. 

Comp. Sci. 34 169-206. 

[6] Duchon P 2000 On the Enumeration and Generation of Generalized Dyck Words. Disc. Math. 
225(1-3) 122-135. 

[7] Dvoretzky A and Motzkin T 1947 A Problem in Arrangements. Duke Math. J. 14 305-313. 
[8] Feller W 1968 An Introduction to Probability and its Applications. Wiley (New York). 
[9] Flory P J 1955 Statistical Thermodynamics of Semi-Flexible Chain Molecules. Proc. Roy. Soc. 
(London) A 234 60-73. 

[10] Flory P J 1969 Statistical Mechanics of Chain Molecules. Wiley Interscience (New York). 
[11] Gessel I M 1986 A Probablistic Method for Lattice Path Enumeration. J. Stat. Plan, and Interfer. 
14 49-58. 

[12] Hammersley J M and Morton K W 1954 Poor Man's Monte Carlo. J. Roy. Stat. Soc. B 16 
23-38. 

[13] Hammersley J M and Whittington S G 1985 Self-avoiding Walks in Wedges. J. Phys. A: Math. 
Gen. 18 101-111. 

[14] Hillc E 1948 Functional Analysis and Semi-Groups. AMS Colloquim Publications 31. American 

Mathematical Society (Providence, Rhode Island). 
[15] Janse van Rensburg E J 1999 Adsorbing Staircase Walks and Polygons, Annals of Combinatorics 

3 451-473. 

[16] Janse van Rensburg E J 2000 The Statistical Mechanics of Interacting Walks, Polygons, 

Animals and Vesicles. Oxford Lecture Series in Mathematics and its Applications 18. Oxford 

University Press (Oxford). 
[17] Janse van Rensburg E J 2005 Square Lattice Directed Paths Adsorbing on the Line Y = qX, 

Journal of Statistical Mechanics: Theory and Experiment (September 2005), P09010. 
[18] Janse van Rensburg E J 2005 Adsorbing Bargraph Paths in a q- Wedge, Journal of Physics 

A: Mathematical and General 38, 8505-8525. Corrigendum (2006): Journal of Physics A: 

Mathematical and General 39, 3847. 
[19] Janse van Rensburg E J 2006 Forces in Motzkin Paths in a Wedge, Journal of Physics A: 

Mathematical and General 39, 1581-1608. 
[20] Janse van Rensburg E J and Le Y 2005 Forces in Square Lattice Directed Paths in a Wedge, 

Journal of Physics A: Mathematical and General 38, 8493-8503. 
[21] Janse van Rensburg E J, Prellberg T and Rechnitzer A 2007 Partially Directed Paths in a 

Wedge. To appear in J. Comb. Theory, Ser. A. 
[22] Madras N and Slade G 1993 The Self- Avoiding Walk. Birkhauser (Boston). 

[23] Mishna M and Rechnitzer A 2007 Two Non-Holonomic Lattice Walks in the Quarter Plane. 

Arxiv Preprint math.CO/0701800. 
[24] Petkovsek M and Bousquet-Melou M 2000 Linear Recurrances with Constant Coefficients: the 

Multivariate Case. Disc. Math. 225(1-3) 51-75. 
[25] Prellberg T and Brak R 1995 Critical Exponents from Non-Linear Functional Equations for 

Partially Directed Cluster Models. J. Stat. Phys. 78 701-30. 
[26] Prellberg T and Owczarek A L 1995 Stacking Models of Vesicles and Compact Clusters. J. Stat. 

Phys. 80 755-79. 

[27] Privman V, Forgacs G and Frisch H L 1988 New Solvable Model of Polymer Chains Adsorbing 

on a Surface. Phys. Rev. B 37 9897-900. 
[28] Sec http://www.research.att.com/ nj as/sequences/ 



